Search for: All records

Abstract Recent marine heatwaves in the Gulf of Alaska have had devastating impacts on species from various trophic levels. Due to climate change, total heat exposure in the upper ocean has become longer, more intense, more frequent, and more likely to happen at the same time as other environmental extremes. The combination of multiple environmental extremes can exacerbate the response of sensitive marine organisms. Our hindcast simulation provides the first indication that more than 20% of the bottom water of the Gulf of Alaska continental shelf was exposed to quadruple heat, positive hydrogen ion concentration [H⁺], negative aragonite saturation state (Ω_arag), and negative oxygen concentration [O₂] compound extreme events during the 2018–2020 marine heat wave. Natural intrusion of deep and acidified water combined with the marine heat wave triggered the first occurrence of these events in 2019. During the 2013–2016 marine heat wave, surface waters were already exposed to widespread marine heat and positive [H⁺] compound extreme events due to the temperature effect on the [H⁺]. We introduce a new Gulf of Alaska Downwelling Index (GOADI) with short‐term predictive skill, which can serve as indicator of past and near‐future positive [H⁺], negative Ω_arag, and negative [O₂] compound extreme events near the shelf seafloor. Our results suggest that the marine heat waves may have not been the sole environmental stressor that led to the observed ecosystem impacts and warrant a closer look at existing in situ inorganic carbon and other environmental data in combination with biological observations and model output. 

Abstract Uptake of anthropogenic carbon dioxide from the atmosphere by the surface ocean is leading to global ocean acidification, but regional variations in ocean circulation and mixing can dampen or accelerate apparent acidification rates. Here we use a regional ocean model simulation for the years 1980 to 2013 and observational data to investigate how ocean fluctuations impact acidification rates in surface waters of the Gulf of Alaska. We find that large-scale atmospheric forcing influenced local winds and upwelling strength, which in turn affected ocean acidification rate. Specifically, variability in local wind stress curl depressed sea surface height in the subpolar gyre over decade-long intervals, which increased upwelling of nitrate- and dissolved inorganic carbon-rich waters and enhanced apparent ocean acidification rates. We define this sea surface height variability as the Northern Gulf of Alaska Oscillation and suggest that it can cause extreme acidification events that are detrimental to ecosystem health and fisheries. 

Abstract. The coastal ecosystem of the Gulf of Alaska (GOA) is especially vulnerable to the effects of ocean acidification and climate change. Detection of these long-term trends requires a good understanding of the system’s natural state. The GOA is a highly dynamic system that exhibits large inorganic carbon variability on subseasonal to interannual timescales. This variability is poorly understood due to the lack of observations in this expansive and remote region. We developed a new model setup for the GOA that couples the three-dimensional Regional Oceanic Model System (ROMS) and the Carbon, Ocean Biogeochemistry and Lower Trophic (COBALT) ecosystem model. To improve our conceptual understanding of the system, we conducted a hindcast simulation from 1980 to 2013. The model was explicitly forced with temporally and spatially varying coastal freshwater discharges from a high-resolution terrestrial hydrological model, thereby affecting salinity, alkalinity, dissolved inorganic carbon, and nutrient concentrations. This represents a substantial improvement over previous GOA modeling attempts. Here, we evaluate the model on seasonal to interannual timescales using the best available inorganic carbon observations. The model was particularly successful in reproducing observed aragonite oversaturation and undersaturation of near-bottom water in May and September, respectively. The largest deficiency in the model is its inability to adequately simulate springtime surface inorganic carbon chemistry, as it overestimates surface dissolved inorganic carbon, which translates into an underestimation of the surface aragonite saturation state at this time. We also use the model to describe the seasonal cycle and drivers of inorganic carbon parameters along the Seward Line transect in under-sampled months. Model output suggests that the majority of the near-bottom water along the Seward Line is seasonally undersaturated with respect to aragonite between June and January, as a result of upwelling and remineralization. Such an extensive period of reoccurring aragonite undersaturation may be harmful to ocean acidification-sensitive organisms. Furthermore, the influence of freshwater not only decreases the aragonite saturation state in coastal surface waters in summer and fall, but it simultaneously decreases the surface partial pressure of carbon dioxide (pCO2), thereby decoupling the aragonite saturation state from pCO2. The full seasonal cycle and geographic extent of the GOA region is under-sampled, and our model results give new and important insights for months of the year and areas that lack in situ inorganic carbon observations. 

Abstract We demonstrate a linking of moderately high resolution (1 km) terrestrial hydrological models to a 3‐D ocean circulation model having similar resolution in the northern Gulf of Alaska, where a distributed line source of freshwater runoff exerts strong influence over the shelf's hydrographic structure and flow dynamics. The model interfacing is accomplished via mass flux boundary conditions through the ocean model coastal wall at all land‐ocean adjoining grid cells. Despite the high runoff volume and lack of a coastal mixing estuary, the implementation maintains numerical stability by prescribing depth invariant and surface‐intensified inflows at fast and slow discharge grid cells, respectively. Based on comparisons against in situ hydrographic data, the coastal sidewall mass flux boundary condition results in more realistic hindcast surface salinity and salinity gradient fields than models that distribute coastal runoff in the form of spatially distributed precipitation. Correlations with observed thermal and haline monthly anomalies reveal statistically significant hindcast temporal variability during the freshet season when the signal‐to‐noise ratio is large. Comparisons of ocean models forced by high‐ and low‐resolution hydrological models reveal differences in salinity, surface elevation, and velocity fields, highlighting the value and importance of accurate coastal runoff fields. The model results improve our understanding of the regional influence of runoff on sea level elevations and the distribution and fate of fresh water. Our approach has potential applications to biogeochemical modeling in regions where distributed line source freshwater coastal discharges deliver heat, momentum, and chemical constituents that may influence the marine carbon pump.